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We investigate quantum electrodynamic effects under the influence of an external, time-dependent electro- 
magnetic field, which mediates dynamic modifications of the radiative corrections. Specifically, we consider 
the quantum electrodynamic vacuum-polarization tensor under the influence of two external background fields: 
a strong laser field and a nuclear Coulomb field. We calculate the charge and current densities induced by a 
nuclear Coulomb field in the presence of a laser field. We find the corresponding induced scalar and vector po- 
tentials. The induced potential, in first-order perturbation theory, leads to a correction to atomic energy levels. 
The external laser field breaks the rotational symmetry of the system. Consequently, the induced charge density 
is not spherically symmetric, and the energy correction therefore leads to a "polarized Lamb shift." In partic- 
ular, the laser generates an additional potential with a quadrupole moment. The corresponding laser-dressed 
vacuum-polarization potential behaves like 1/r 3 at large distances, unlike the Uehling potential that vanishes 
exponentially for large r. The energy corrections are of the same order-of-magnitude for hydrogenic levels, ir- 
respective of the angular momentum quantum number. The induced current leads to a transition dipole moment 
which oscillates at the second harmonic of the laser frequency and is mediated by second-order harmonic gen- 
eration in the vacuum-polarization loop. In the far field, at distances r ^> 1/lo from the nucleus (lu is the laser 
frequency), the laser induces mutually perpendicular electric and magnetic fields, which give rise to an energy 
flux that corresponds to photon fusion leading to the generation of real photons, again at the second harmonic 
of the laser. Our investigation might be useful for other situations where quantum field theoretic phenomena are 
subjected to external fields of a rather involved structure. 

PACS numbers: 12.20.Ds,11.10.Gh,03.70.+k 



I. INTRODUCTION 

This article is concerned with the calculation of the 
vacuum-polarization tensor in the presence of two nontrivial 
background fields: a nuclear Coulomb field and a (possibly 
strong) laser field. Particular emphasis is laid on the induced 
charge density and potential, as well as the induced current 
and the vector potential, and the resulting induced electromag- 
netic fields. As an application, we consider the energy shifts 
of atomic states in hydrogenlike systems, and induced transi- 
tion probabilities, under the influence of the laser-dressed vac- 
uum polarization. However, we would like to emphasize that 
the two latter aspects constitute merely example applications 
of the more general expressions derived in this article. 

Before we indulge into this endeavour, let us briefly re- 
call a few very well known facts about vacuum polarization. 
It is well known that the Coulomb law becomes invalid at 
distances comparable to the electron Compton wavelength 
A c = 3.8 x 10 -11 cm. This phenomenon is due to the excita- 
tion of virtual electron-positron pairs from the vacuum, which 
polarize the vacuum in a way approximately equivalent to a 
dielectric medium. Indeed, at large distance, we see only a 
"screened" charge, which is less than the "bare" charge visible 
to an observer who approaches the particle in question, in our 
case mostly an atomic nucleus, to a very close distance. The 
total charge density induced in the vacuum has to be zero, and 
this also defines the renormalization condition for the charge 
which can be defined via the low-energy limit of the Thom- 
son scattering cross section (see 1 1], Chap. 7), as it affects two 
distant charge particles. However, as we enter the cloud of the 



electron-positron pairs that surround the charged nucleus, the 
effective charge becomes larger, because one has to add to the 
"renormalized", physical charge of the nucleus (seen at large 
distance), a contribution due to the electron-positron pairs that 
surround the nucleus at short distance. 

So, the virtual excitations of the electron-positron (or 
muon-antimuon) pairs naturally generate a "fifth-force" like 
Yukawa potential contribution to the Coulomb law, which 
is however exponentially suppressed at distances larger than 
the Compton wavelength, and by a quantum electrodynamic 
(QED) perturbation theory parameter a. 

What happens to the vacuum-polarization effects in a laser 
field? Evidently, the laser field, which is composed of many 
photons, will modify the electron-positron loop as the laser 
photons interact with the virtual quanta during their short life- 
time that is only of the order of X c /c = 1.3 x 10 -21 s. We find 
that in a remarkably strong laser field, the modification of the 
vacuum-polarization contribution to the Coulomb law can be 
significant. In particular, the laser induces a rotational symme- 
try breaking in the correction to the Coulomb law, which leads 
to a polarized Lamb shift whereby states with different mag- 
netic quantum numbers receive different energy shifts. The 
induced charge density, in a laser field, has a monopole and a 
quadrupole component. The quadrupole part of the potential 
is not exponentially suppressed at large distances, but tails off 
like a "normal" quadrupole potential, i.e. as 1/r 3 . 

In part, our investigation follows ideas outlined previously 
in the context of the analysis of the vacuum-polarization ten- 
sor in a laser field Here, our focus is on applications 
including a (possibly strong) binding Coulomb field in addi- 
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tion to the laser field. Another fundamental QED process, the 
self-energy under the influence of an external laser field, for a 
free particle, has previously been studied in 0,0] .and for a 
strongly driven transition between bound states, in • 
A more general context of our study might be defined as 
"quantum field theory under the influence of external fields." 

From a naive point of view, one might expect that in a 
strong laser field, the charge density induced by the vacuum 
polarization loses the rotational symmetry, and that an align- 
ment along the laser propagation and polarization axes oc- 
curs. However, the exact functional form of the induced laser- 
dressed contribution to the charge density [see Eq. J26l > be- 
low] is not accessible by symmetry arguments. It is therefore 
indispensable to use a fully laser-dressed formalism, where 
the quantum electrodynamic effects are incorporated into the 
propagators right from the start of the calculation. As will be 
discussed below, the result of a fully laser-dressed calculation 
confirms the intuitive picture, according to which rotational 
invariance is broken in both linear as well as circularly polar- 
ized laser fields, and a few unexpected asymptotic structures 
are found for the laser-induced fields. 

This article is organized as follows: in Sec. |H] we analyze 
the additional charge density induced by the external laser 
field (this potential is an addition to the "normal" Uehling 
potential). The way from the charge density to the induced 
potential is detailed in Sec. Ill Al The asymptotics of the in- 
duced potential, in various regions from the atomic nucleus, 
are derived in Sec. Ill Bl As an application, the energy shift of 
hydrogenic bound states is considered in Sec. Ill CI 



II. INDUCED CHARGE DENSITY AND POTENTIAL 

A. From the charge density to the potential 

The amplitude T of photon scattering in the laser field has 
the form 

T = 4*T fiV (k 2> k 1 )e v 1 . (1) 

Here, the four-vectors fci and k 2 characterize the incoming 
and outgoing photon momenta, as they flow into and out of 
the laser-dressed loop. The polarization vectors of the pho- 
tons are e^ 1 and e^. The quantity T IJi ^(k 2 , fci) is canonically 
identified as the vacuum-polarization tensor. The calculation 
of the tensor T )JiU (k 2 , fci) is a rather involved calculation, but 
complete results, valid for arbitrary laser-field configurations, 
were found with the use of operator techniques in Ref. Q|. 
Independently, by means of another method, the vacuum- 
polarization tensor was obtained in Ref. |3]. Both results are 
in agreement with each other. Of course, the scattering, which 
is beyond Maxwell's classical electrodynamics, appears due 
to the polarization of the electron-positron vacuum in the pres- 
ence of a laser field. 

In this article, we use rationalized Gaussian units with h = 
c = 1, and eo = (47r) -1 , which implies e 2 = a where e is 
the electron charge and a is the fine-structure constant. We 
represent the vector potential of the laser field in the form 










FIG. 1: Diagrams corresponding to the amplitude < I 1 1 . including 
their multiplicities. The solid lines are the electron propagators, the 
dashed line with the cross denotes the Coulomb field, the wavy lines 
correspond to the laser photons, and the zigzag line is the induced 
field. 



Here, <f> = xx = x°t — xfis the phase of the laser field at 
the space-time point (r,t), and u = xP = \x\ is the laser 
frequency. We use the spacetime metric in the Bjorken-Drell 
convention g M „ = diag(l, — 1, —1, —1). The vector x points 
into the laser propagation direction, and ipi^irf 1 ) are functions 
which characterize the shape of the laser pulse. We denote the 
(orthogonal) unit vectors giving the polarization directions by 
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ei and e 2 , respectively (ef = e 
four- vectors e\ = (0, e*i) and e 2 



— 1). The corresponding 
(0, e 2 ) satisfy 



x 2 = x ■ e\ = x ■ e 2 = e\ ■ e 2 = , e\ = e\ = — 1 . (3) 

All vector scalar products (involving three or four- vectors) are 
denoted by • in this article. We also use the conventions 



v = — , n = — 

to r 



(4) 



for the unit vector v — e i x e 2 pointing into the propaga- 
tion direction, and the unit vector pointing to a coordinate r, 
respectively. 

In Refs. 12, y|], the photon scattering amplitude was found 
for arbitrary functions ?/>i,2- In the current paper, we restrict 
ourselves to the case of monochromatic plane wave with the 
frequency u>, although with allowance for arbitrary (linear, 
elliptic) polarizations. In this case, the functions tpi((f>) and 
ip 2 (4>) take the form 



ipi (</>)= a\ cos <ft , ip 2 (<j)) — a 2 sin <p , x = t 
The scattering amplitude depends on the parameters 
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where m is the electron mass. The mean square laser field 
strength £ 2 , averaged over a laser period, is expressed via the 
parameters Q 2 by the relation 



= ei M Vi(0) + e 2 ^ 2 {(j)) . 



(2) 
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FIG. 2: (Color online.) The familiar rotationally symmetric vacuum- 
polarization charge density around an atomic nucleus is severely dis- 
torted in a strong laser field. Here, we plot the additional laser- 
dressed, time-independent contribution to the renormalized charge 
density Sp R (r), as induced by the laser field. The quantity 8p^(r) 
has to be added to the well-known expressions for the Uehling and 
the Wichmann-Rroll parts (for a review see e.g. fl(il ). The result for 
Spo (f) is obtained after Fourier transformation of the momentum- 
space result given in Eq. d 1 8i below, into coordinate space. Along the 
vertical axis, we have the propagation direction of the laser, which we 
assume to be linearly polarized (£1 7^ 0, £2 = 0). The polarization 
direction (laser electric field) is along the a;-axis. We plot a charac- 
teristic surface given by a constant induced charge density, namely 



the surface defined by the relation Spo(r) — |2 



(S/SoY K, 



where K is a numerical constant given by K = ±1.204 for the 
plot shown. Two areas with a positive induced charge density (red, 
aligned mainly along the y axis, K = +1.204) and a negative one 
are distinguished (the latter is in blue color and located along the 
propagation axis, K = —1.204). As emphasized in the text, the total 
induced charge is zero. 



L 



L 




(b) 




(c) 



Thus, we may say that the parameters £x,2 measure the laser 
field strength. The parameter A, by contrast, measures the mo- 
mentum scales relevant to the vacuum-polarization loop (in- 
coming momenta and the four-momentum of the laser wave) 
against the electron mass scale. 

The following identity is sometimes useful: 




Here, £q is the Schwinger critical electric field strength £rj = 
m 2 /e = 4.4xl0 15 Vcm" 1 . 

As is well known, the presence of an external field leads 
to the existence of an induced current J p due to the vacuum 
polarization. If the external field is a pure Coulomb field, 
then for the induced current we have J = 0, and p = J° 
is well known in this case (see, e.g., lUlll ). Here, we consider 
the situation where the external field is a Coulomb field, but 
there is, in addition, the laser field, which dresses the vacuum- 
polarization loop. 

Our incoming photon is the Coulomb quantum with ki = 
q = (0,q) and = 6g. Then, J* = T> l0 (k 2 ,k 1 ). We 
calculate J M in the leading approximation in the parameter 
Za, where Z is the nuclear charge number. Due to momentum 



FIG. 3: (Color online.) Projections of the laser-induced part of the 
vacuum-polarization charge density looked upon from the x, y and z 
axis, in Figs, (a), (b) and (c), respectively. The conventions are the 
same as in Fig. [2] the z-axis is the laser propagation direction, and 
the x-axis is the polarization axis (laser electric field). The picture 
provides a ray-traced impression of the charge density, represented as 
an absorptive medium. Red color is used for positive induced charge 
(areas located along the x and y axes), blue color is used for negative 
charge (areas located along the z axis). The plots are as seen by an 
observer looking at the atomic nucleus from the specified directions. 
This is in contrast to Fig. |2| where the isosurface corresponding to 
a given, constant charge density is plotted. In the above graphical 
representations, the fall-off of the induced charge density for large 
radial arguments can be distinguished clearly. 



and C-parity conservation (Furry theorem), we have k 2 — 
k% + 2lx, where I — 0, ±1,±2..., and this fact affords an 
explanation for the two alternative definitions of A in Eq. 0. 
A noticeable contribution to the induced charge density arises 
at distances r ~ A c = X/m , which reciprocally corresponds 
to the momentum of a Coulomb quantum \q\ ~ m. 

For all current high-intensity lasers, the laser photon energy 
is small as compared to the electron rest mass energy (ui <C 
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m), and therefore 



|A| 



\q\u> u> 

< — < 1 



(9) 



Our calculations as reported in the current article depend on 
the approximation A 1 [see also Eqs. (2.5), (2.31) and 
(2.32)ofRef. |2]]. All results reported below are given on the 
first nonvanishing order in the expansion in powers of oj/m. 
A further condition imposed on our calculations concerns the 
laser intensity parameters £i i2 , 



< e h2 « 1/A 2 

which according to Eq. (|8j means that 



to 
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The asymptotics of T 00 (fc 2 , fa) = T 00 (fc 2 , q) under this ap- 
proximation are found to consist of terms with absorption and 
emission of one laser photon each, and of terms with absorp- 
tion of two and emission of two laser photons (see Fig.^, 



T 00 (fc 2 , g ) = (2^) 4 [S(q-k 2 )B^(q) 

+6(q - k 2 - 2x) (q) + 5(q - k 2 + 2x) (q) 

(11) 



in the matrix element. Here, the functions B^\q) and 
(q) are given by 



B<-°Hq) = 



aUJ " ' ' '- ,.2\2 n 2 J 3. 2 t2\/„2 „2\ 



(e 2 +£ 2 2 ) 

(1-u 2 ) 



( g 2 + ^(l + , 2 + £)+ig(l- w 2 ) 2 £ 



(12a) 



1 

B{±){<f) = 4^ /^(i-- 2 ) 2 ^ 2 {±3 i 66gig 2 + ^« 2 +e 2 2 )(g 2 - 92 2 ) 



(1-v 2 ) 



(9 



2 \2 



(rf +(1 g)(l+«»+2J) + ^-(l_^)»I> 

4m z 



(12b) 



where 



<7i = q ■ ei, g 2 = 9 



£> = £(q 2 ) = 



1 



e2, 
' 2 (1 
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In the derivation of Eq. d!2l >. we have neglected higher-order 
terms in the expansion in powers ofui/m. The Fourier trans- 
form of the correction to the induced charge density due to the 
laser field, 6p(q, t), is essentially given by the Fourier trans- 
form of T 00 (fc 2 , q). It is a sum of time-dependent and time- 
independent parts: 



Sp(q,t) = Sp (q) + Spx(q,t), 



(14) 



where 



Ze 
Ze 



is (0) (9), 



(15) 



^-\q) 



(q- 22) 



(«f+2x) 2 ' 

We now turn to a more specific evaluation of the charge den- 
sity, in various regions of practical interest. First, we investi- 
gate (in coordinate space) the distance r ~ ao = l/(Zam), 



which typically gives rise to the main corrections to the en- 
ergy in bound systems (do is the Bohr radius). For a laser 
whose frequency is of the order of atomic transitions, we have 
uj ~ (Za) 2 m <C 1/clq, so that the "interesting" distance ful- 
fill r <C 1/w, and the corresponding momenta are in the range 
q ~ 1/r 3> (J = We can therefore replace 



(q± 2k) 2 



(16) 



in Eq. dl5l . Note that this replacement, in particular, does not 
affect the validity of the expressions for radial arguments of 
the order of the Compton wavelength. Employing the approx- 
imation ( I16> . we obtain from Eqs. (I12> . (I14> and d!5> 



rV/> 2 



247T m 2 q 



4 

if' 2 
4m 



(qi-qim-e 2 ) + 



{ -P^(1 + v 2 + D) 
(1 — v l ) 



+^zf(i-« 2 ) £) 



(17a) 
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and 



a(Ze)uj 2 
247T m 2 q 2 



dv(l-v 2 ) 2 D 2 



(17b) 



Hi - £2) cos(2wi) + 36 qi Q2 sin(2wi) 



Although the expression for the vacuum polarization tensor 
in Eqs. (2.5), (2.31) and (2.32) of Ref. [2] is in principle al- 
ready renormalized, it is still necessary to carry out an addi- 
tional finite subtraction, namely, to subtract the asymptotics 
of the induced charge at \q\ — ► 0. This is akin to the ad- 
ditional subtractions that have to be carried out in the calcu- 
lation of the Wichmann-Kroll correction. The renormaliza- 
tion condition in the on-mass-shell scheme, for the Coulomb 
field is that two distant observers should observe exactly the 
Thomson scattering cross section, which implies that the total 
induced charge has to be equal to zero (an eludicating discus- 
sion surrounds Eq. (7.18) of Ref. 1 1]). In addition, our sub- 



traction provides a finite quadrupole moment, corresponding 
to the induced charge distribution. Because 5p\(q,t) — > 
at |<t*| — ^ [see Eq. (ll5H . it is only the time-independent part 
Spo(q) which receives an additional, finite renormalization. 

The subtraction leads to the renormalized time-independent 
component of the charge density Sp R (q), 



s Po(q) = § po(q) 



a(Ze)co 2 
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60-7T m 2 q 2 



e 2 M 



(18) 



The renormalized charge density 8p R (f), excluding the Dirac 
delta contribution, is represented graphically in Figs. [2] 
and [3] The full renormalized Sp R (q,t) (including the time- 
dependent term) is 



5p R (q,t) = 5p R (q) + S Pl (q,t). 



(19) 



The potential $(r, t) in the Lorentz gauge satisfies the equa- 
tion 

|L$(r,t)- V 2 $(r,i) = 4nSp R (f,t). (20) 



Therefore, the Fourier transform of the potential, which corresponds to the density Sp R (q, t), is 



%5p R (q) 



4tt 



Spi(q,t) 



4-7T 

-5p R (q)+^(Ze) 
q- 



Jliujt 



B^(q) 



-2iuit 



(r 



Auj 2 



i0)(q + 2>c) 2 (q 2 -4oj 2 -i0)(q-2x) 2 



(21) 



Here the regularization of the denominators, (±i0), corresponds to the boundary condition of radiation (see below). It is conve- 
nient to present the potential in coordinate space as a sum of time-dependent and time-independent parts: 



$(r,i) = $ (r) + $i(r,i). 



(22) 



By Fourier transform of Eq. i2l\ . we thus obtain at a distance r <C 1 /w from the nucleus, 



Mr) = 



$!(r,t) 



a(Ze)uj 2 



[F x {mr) + F 2 (mr) (l - 3(z7- n) 2 )] (g + 



487rm 

-F s (mr) [{fi-e x f-{n-e 2 f] 
a(Ze)uj 2 



487TTO 



Fi(mr) + 



45mr 



F 2 (mr) 



22 
45mr 



'1 - i{vn) 



(23a) 



(tf -$)cos(2wt) 



- \F 3 (mr) - JL} [{£ + e|) cos(2urf) [(ft ■ e\) 2 - (ft ■ e 2 ) 2 ] + 4£i6 sin(2a,i)(n ■ e x )(ft ■ e 2 )] j, (23b) 
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where we note the conventions in Eqs. (|3j and (0}. In this formula, 

1 



Fi (mr) 
F2(mr) 



dv (1-v 2 ) 1 / 2 1.^(1 -v 2 )(l + 2y) - 
dv(l-v 2 )^ 2 (^- 



4 

3y 



(l + v 2 )(l + y) + l + ^y+^y 2 



-2y 



(24a) 



(l + v 2 )(l + y) + l + ^y+^y 2 



-2y 



+ 1 -±^[l-(l + y )(l + y + y 2)e- 2 y]+ 



F 3 (mr) 



9 f , (1-w 2 ) 3 / 2 
dv - 



y = y{r, v) = 



y 

mr 



2y3 
2 



1 - (1 + 2y + 2y 2 + 7 -y 3 + \y 4 )e~ 2 y 
6 3 



(24b) 



l-(l + y)(l + y+^y 2 )e- 2y 



s/T 



(24c) 



r 



In principle, the expressions i23\ gives the complete answer 
for the induced potential, which is necessary for the calcula- 
tion of the energy shift. We reemphasize that the only addi- 
tional approximation employed was A 1, which is easily 
fulfilled for current and planned high-intensity laser facilities. 
Note that in the particular case of the circular polarized laser 
field (£i = ±£2 = 0' tri e expression for the induced potential 
becomes essentially simpler. 



B. Asymptotics of the induced potential 

Let us consider the behaviour of the potential d23t at large 
and small distances. The functions Fi(p) have the following 
asymptotics at p = mr 3> 1: 



1/m. We recall that the representation in Eqs. d23l and &241 
is valid under the (additional) condition r <C 1 /u>. Thus, we 
write the asymptotics of the potential atl/cj^>r^> 1/m in 
the form 



<Z>(r,t) 
f(n,t) 



rtniQV 1 



(26) 



a(Ze)u) 2 



llh + ^-n) 2 ] (£ 2 -£ 2 2 )cos(2^) 



3607T m 2 

+ml + CI) cos(2^) [(ft ■ el) 2 - (n ■ e 2 ) 2 ] 
+12£i£ 2 sin(2wi) (n • e\) (n ■ e 2 ) 



Flip) 



(i 



e~ 2p , F 2 (p) 



148 
105p 3 



F 3 (p) 



36 
35p 



•(25) 



We can thus safely neglect the exponentially decaying Fi (p)- 
term against the polynomial behaviour of F2 and F3 for r ^S> 



(27) 



The first term on the right-hand side of Eq. i26\ . proportional 
to 1/r 3 , can be associated in a natural way with a charge dis- 
tribution that originates from an oscillatory quadrupole mo- 
ment, 



'Sa(Ze)u)' 

357TTO 4 



- jg (6 ij - 3I/V) sin 2 ^) + & cos 2 ^)] 



-{e\e{ - e 2 e 2 ) [£ 2 sin 2 (wt) - ff cos 2 (ut)] - (e\e J 2 + e 2 e{)£i£ 2 sin(2wi) 



(28) 



The second term on the right-hand side of i26\ . proportional 
to 1/r, reflects the existence of radiation in the system, but 
cannot be used for a valid description of the radiation. The 
reason is that the derivation of the asymptotics (I26> was car- 
ried out at distances 1/m <C r <C 1/w. This is insufficient for 
the consideration of radiation, which by definition is only a 



well-defined process at distances much larger than the wave- 
length of the emitted light. For the calculation of the radiation, 
it is necessary to consider distances r ^> 1/ui (see Sec. lllIDl 
below. 

As was already pointed out, for the purpose of calculating 
the energy shift, it is necessary to know the time-independent 
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part of the potential (for the oscillatory terms, the time- 
averaging is applied over a full period of the driving laser). 
The time-averaged part of the quadrupole moment has the 
form: 



707TTO 4 [27 



(29) 



At small distances from the nucleus, r <C 1/m, we get 
a complete different asymptotic behaviour for the functions 

Fi{p), 

oq 99 O 

Fl („),-_, M*—, *(,)*-. (30, 
For the potential, this implies 



3607rm 2 r 

+3[(n-e 1 ) 2 ~(n-e 2 ) 2 } - £ 2 2 ) \- (31) 



We see that the leading small-distance asymptotics is indepen- 
dent of time. 



C. Application: energy shift of hydrogenic levels 

Let us consider the correction 8E = (e&o(r) ) to hydro- 
genic energy levels, corresponding to the potential $o(r), 
Eq.(l23>. For n s^-states, only the spherically symmetric part 
of this potential gives a nonzero contribution. It is propor- 
tional to the function Fi(mr). Because the function Fi(mr) 
decreases exponentially at mr ^> 1, we can replace in the ma- 
trix element the wave function ip(r) by if>(0), because it has 
the electron Compton wavelength is much smaller than the 
characteristic scale ao = l/(Zam) 3> 1/m. 

With the considerations, the energy shift for S states can be 
easily calculated, and we obtain [see Eqs. (|8j, (I23i and (I25H 



SE(nS) 



8a(Za) 4 cj 2 
1357rn 3 m 

16a(Za) 4 m ( £ 

1357rn 3 \£q 



(e 2 + e 2 ) 



(32) 



Therefore, due to the factor £/£q, in all realistic cases the 
correction obtained is small. 

Let us consider the corrections to the energy level with the 
orbital angular momentum I > 0. In the general case of an 
elliptically polarized wave (£i ^ £2) , the projection j z of the 
total angular momentum on the propagation direction of the 
laser v is not conserved. This implies a nonvanishing nondiag- 
onal matrix element of the potential <f>o (r) for Aj z — ±2. For 
j > 3/2 (where j is the total electron angular momentum), 
it is necessary to solve a more complicated secular equation 
in order to find the eigenvalues in the presence of the laser- 
dressed vacuum polarization. For simplicity, we restrict our- 
selves to the case of a circular polarized laser field, where the 



mentioned problem does not appear, and all matrix elements 
remain diagonal. 

For I 0, the main contribution to the integral for the 
matrix element &o(r) comes from distances of the order of 
the Bohr radius, r ~ ao 3> 1 /m, and is generated by the 
quadrupole potential. Therefore, we can substitute to the 
matrix element the asymptotics of the potential < f ) o(r) ~ 
Qq 3 n % v? /(2r 3 ) with the quadrupole moment [see Eqs. J26l > 
and ( I29H being equal to 



37a(Ze)uj 2 Z 2 
315?r m 4 



{5 l - 



3i/V 



(33) 



for the case of circular polarization of the laser field. Using 
the relations fTHl 



{Zaun) 



'•1(1 



dn\Y lm \ 2 (i 



-1)(H 
3 cos 2 



1/2)' 



3m 2 - 1(1 + 1) 
(2Z + 3)(/-l/2) ' 



(34) 



where I is the orbital angular momentum, we obtain the fol- 
lowing shift of energy level with I ^ 



6Ej 



yia(Za) 1 m 
3157m 3 



3j?-j(j + l) 
1(1 + l)(2l + l) 2 



°j-i/2,l 



J + l 



j 



(35) 



In rewriting the result in terms of the total rather than orbital 
angular momentum of the electron, we have used the explicit 
representation of the spinor spherical biharmonics II ill . It is 
interesting that the magnitude of the correction for I ^ is 
of the same order as that for I = 0. This is in sharp con- 
trast to the "normal" vacuum polarization in a Coulomb field, 
where the dominant effect is given by a potential of the form 
— (1/15) aZam 2 5(r), which has a nonvanishing expecta- 
tion value only for S states. 



HI. INDUCED CURRENT AND VECTOR POTENTIAL 
A. From the current to the vector potential 

In a pure Coulomb field acting as an external "condition", 
or source, it is well known that the spatial components of 
the induced current vanish completely. The situation changes 
drastically for an additional laser field. Then, the spatial com- 
ponents of the induced current no longer vanish, and that 
constitutes an additional effect supplementing the vacuum- 
induced laser-dressed charge density Sp, which exists in both 
cases. We discuss the Fourier transform of the induced non- 
vanishing current density J(q, t). We employ analogous steps 
as in the derivation of Eq. (II \\ . 

We define the following vector to represent the mixed 
spatial-temporal components of the vacuum-polarization ten- 
sor, 



T^(k 2 ,k 1 ) = T\k 2 ,k 1 ). 



(36) 
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In analogy to Eq. i2l\ . we find the following structure, 



The functions G {0) (q) and (q) are found to read 



f(k 2l k 1 ) = {2nf[8{k 1 -k 2 )& a \q) 

+ 8(ki - k 2 - 2x)G { -\q) + 5{ki -k 2 + 2x)G i+) {q)} 



(37) 



G (0) (g) 



24-7T m? 



487r m 2 



dv (1 - v 2 fD 2 \ f±f (2 + v 2 D) [qx \V x q]} - h& - g) [qx [V x {q x 



92 J 





^M.( 2 + v 2 D) [±2w(# 



dv (1 - </)^D <| -(£ + C|) [±2w(& - g 2 ) + [?x p x (ft 



92 J 



[g x [v x q]]] ± - i £i£ 2 [{±2ujv -q)x(q\- q 2 )} [V ■ (ei x e 2 )] 



Here, D is defined as in Eq. Jl 3i . and we recall the conditions 
l|9} and dl 01 which apply to all results reported in this paper. 

Although v = e i x e* 2 an d ? 2 = 1 in a specified coordinate 
system [see Eq. (0)], it is still useful to indicate the factor 



v- (ei x e 2 ) 



(39) 



The reason is that under a parity transformation, G^\q) 
should transform like a vector, and the expression (±2wi? — 
q) x (qx—qV) would otherwise transform like an axial vector. 
In Eq. J38I . D is defined as in Eq. Jl 31 . and q\ = ei{eiq). The 
result is valid under the conditions l|9} and dlOi . 

In analogy to Eq. dl4> . we represent the induced current as 
a sum of the time-independent and time-dependent parts: 



where 



Mo) 

Ze 
Ze 



J(q,t) = J (q) + Ji(q,t) , 
&°Hq) 



(40) 



(41) 



(q + 2k)' 



:G { -\q) 



-liwt 



(q- 2k)' 



-G (+) (q) 



In analogy to Eq. Jl 81 . we now have to subtract from Jo (q) 
its asymptotics at \q\ — > 0. A time-independent component of 
Jo(q), constant in momentum space and therefore local in co- 
ordinate space, would otherwise lead to a magnetic-monopole 
field. In order to preserve the nonexistence of a magnetic 

monopole, we introduce a finite subtraction J(q) in the time- 
like component of the current Jo(q), so that the renormalized 
current J R (q, t) has the form 

J*(g, t) = Jj(q) + J x {q, t) = [J (q) - J(q)} + J t (q, t) . 

(42) 



(38) 



I 

Here, the finite renormalization J(q) has the form 
a(Ze)co 2 



1807rm 2 <f 2 
-3(Ci 2 -Cl)[3*x [vx(q 



(43) 



In order to infer the induced vector potential A(f, t), one 
now has to solve the equation 



0_ 

dt 2 



A(r, t) - W 2 A(r, t) = 4tt J R (r, t). 



(44) 



In Fourier space, the solution is immediate, and A(q, t) is 
found to be 



A(q,t) =%j£{q) 



4ir(Ze) 



e 2iut G^(q) 



{q 2 - Alj 2 + i0)(q '+ 2k) 2 
2wt &+\q) 



{q 2 - Auj 2 - i0)(q- 2k) 2 



(45) 



B. Current and vector potential in various asymptotic regions 

We first consider the momentum range \q\ ^> u>, which is 
relevant for atomic physics phenomena. In coordinate space, 
this condition corresponds to distances r <C 1/uj. We find 
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l 



o 

-[6 2 sin 2 (^) -£ 2 cos 2 (wt)] [gx [i?x (ft - g 2 )]] - 66 sin(2urf) [?x (gl - g 2 )] (ei x e 2 )]| - J(g) , (46) 

where we recall the definition of D in Eq. (II 31 . The renormalized current J R (q, t) contains both time-independent components 
as well as time-dependent components (the latter vanish when averaged over a full laser period). 

We now use the additional approximation r ^> 1/m. We find that the asymptotics of the induced current density in the range 
1/u) ^> r ^> 1/m are given by the purely time-dependent component 

J*m - Ur,t) - g^f^ {y^i - ® 3 ^ -**)]«»(2art) 

-(fi+fi) -«a) - fa -n 2 )(i7-n)] cos(2urf) + 266 sin(2urf) [n X (rii - n 2 )] (e\ x e 2 )]} . (47) 

Here, we have used the conventions 

rii = (et -n) , Hi — (ei ■ n)e*j . (48) 

Thus, Ji(r, i) cx 1/r 3 in the range 1/uo ^> r ^> 1/m. The time-independent current Jo*(r) at distances r ^> 1/m diminishes 
exponentially, in a way similar to the time-independent induced charge density p R (r), and therefore is omitted in Eq. J47i . 

Again, in the range 1/m < r < we now calculate the time-independent and time-dependent components of the vector 
potential, Ao(r) and A%(f, t), respectively. These correspond to the currents Jq(t) and Ji(r, t). The result reads 

Mr) = Hof^ jf + P- 3n(^ ■ n)] - 3(£ 2 - £ 2 ) [*?(n 2 - n 2 ) - (ft - n 2 )(*7- ft)] ] , (49a) 



and 



a(Ze)uj 



Al (r, t )=~ q fln 2 1 " fi)] cos(2^) + 3(6 2 + f|) [*(n? - n 2 ) - (ft - ft)(i7- ft] cos(2 w i) 

3oU7r m/r ( L 

- 666 x (ni - n 2 )] [v- (S x x e 2 )] sin(2wi) 1 . (49b) 



We emphasize that the behaviour of Ai(r, t) corresponds to the existence of radiation in the system, in analogy to the corre- 
sponding time-dependent term in Eq. (1261 : however, the range of radial arguments 1/m <C r <C 1/w is inappropriate for a 
discussion of this phenomenon. The radial distances relevant for radiation will be discussed below in Sec. lIIIDl 

Using the formula H = V x A, applied to the results in Eqs. (|39a) and ( !49bL it is then very easy to obtain the corresponding 
asymptotics of the magnetic field in the relevant range of distances, 



^o(r) - tn^l^ ~ Snnmi [P ■ (Si x e 2 )] [(fix -n 2 )xv][, (50a) 



+ (£.!+ £,2) nn in2 (ei X e 2 )] cos(2wt) -66" (ei X e 2 )] (rij - ri z 2 ) sin(2urf)j . (50b) 

At small distances, r <C l/m, the leading term of the asymptotics of the vector potential and the magnetic field is time- 
independent: 

^ = ~ ^n e) i \ U & + ® W ■ ")1 + 3 (£ - ^ [r? x \V x (ft - ft)]]| , (51) 
3607rm/r [ J 

a(Ze)oj 2 
1807T m 2 r : 



^o(r) - i a l Ze)tJ ?\ \ U«i + fa) [« x ^ + 9(^ 2 - ff) ffnma [* ■ (Si x e 2 )] ^ . (52) 



I 

In passing, we note that Ho(f) of course has the transforma- necessary introduction of the factor v ■ (ei x e 2 ) in Eq. ( 15 11 1 
tion properties of a pseudovector, which is consistent with the 
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Also note that at all distances, Ho (r) is an odd function of the 
vector ft. Therefore, due to the parity conservation, for any 
hydrogenic atomic state, the energy shift related to the mag- 
netic field vanishes in first order of perturbation theory. 



C. Application: transition probability in atoms 

As shown in previous sections of this article, the laser- 
dressed vacuum polarization induces additional electric and 
magnetic fields which modify the atomic transition probabili- 
ties. We perform the calculation to the first order of pertur- 
bation theory with respect to these fields. The most inter- 
esting situation occurs when the laser frequency u> is closed 
to AE/2, where AE is the energy interval between levels, 
because the second-harmonic generation from the vacuum in 
this case leads to the generation of resonance photons. For 
simplicity, we restrict ourselves by the case of circular po- 
larization and transition Is — > 2p (which corresponds to the 
electric dipole transition for spontaneous radiation). We direct 
the quantization axis along the vector v. Due to the conserva- 
tion of angular momentum, the matrix element is nonzero if 
Aj z = ±2 only. The sign in this selection rule depends on the 
helicity of the laser wave. 

Taking into account the selection rule, let us consider, for 
instance, the case of a right helicity of the laser wave and the 
transition from the state lsx/a with j z = —1/2 to the state 
2?>3/2 with j z — 3/2. Because the induced potential < &i(r, t) 
is an even function of the vector fi, see Eq. d23i . this poten- 
tial does not contribute to the matrix element to the first order 
of perturbation theory, and the leading nonzero contribution 
comes from the perturbation Hamiltonian 



V 



-(Ai(?,t)0) - fi(aHi(r,t)) 



(53) 



where fi = e/(2m) = — \e\/(2m) is the electron magnetic 
moment, and p is the momentum operator. For the vector po- 



J 



2^3/2, jz = 3/2 



-(A(r,t)p) — 
m h. 



tential and the magnetic field, we should use its asymptotics 
at 1/m <r< 1/oj, as given in Eqs. ( I49bl > and J50bi . For s 
states, we can use pipi s (r) = im(Za)nipi s (r), and for cir- 
cular polarization with £ 2 = £f , we have nA\ (r, t) = in the 
asymptotics d49b^ In this setting, only the second term in the 
interaction Hamiltonian in Eq. (I53> gives a nonzero contribu- 
tion. For the latter, a straightforward calculation of the matrix 
element gives 



dt(2p 3/2 ,j z =3/2\ 
a{Zafuj 2 e 



-[1(7 Hi 



|l*i/a,J« = -V2 



6757T771 



2n5{E 2p - E ls - 2uj) , 



(54) 



where Hi = Hi (r , t) and £ 2 = £ 2 = £f f° r circular polar- 
ization. Note that we here assume the quantization axis for 
the total electron angular momentum j to be aligned along 
the laser propagation (not polarization) direction, and the an- 
gular momentum projection j z to be defined with respect to 
that axis. The resonance frequency for a hydrogenic transi- 
tion involving a change in the principal quantum number is 
of order (Z a) 2 , so that the prefactor in Eq. ( I54l > is in fact of 
order a (Z a) 7 £ 2 , where in turn £ 2 is proportional to the laser 
intensity according to Eq. 0. 

The transition Is — > 2p 3 / 2 is driven at the second harmonic 
of the incident laser frequency [see the Dirac delta function 
S(E2 P — Eis — 2a;)], anticipating the possibility of photon fu- 
sion in a strong laser field, to be discussed below in Sec. lIIIDl 
in the far field r> 1/w from the atomic nucleus. 

The transition amplitude in Eq. i54\ should now be com- 
pared to other effects which also can induce small nonvan- 
ishing transition probabilities for lsj/2 — * 2p 3 / 2 , even if the 
incident radiation is off resonance and amounts to only half 
the frequency of the resonance transition. One such term is 
given by 



H 



1*1/2, j. = -1/2 



(55) 



Note that the fields A(r, t) and H(r, t) are the full laser 
fields, i.e. they do not carry any subscripts like the vacuum- 
polarization induced fields. For a transition involving angular 
momentum projections j z = 1/2 — > j z = 3/2 (in contrast to 
j z = —1/2 — > j z —3/2), the term corresponding to d55l with 
two interactions — — (A(f, t)p) may contribute, when we ad- 
ditionally consider higher-order multipoles of the laser field. 
Note that the laser field uniquely specifies a polarization axis 
along which the exponential cxp(ikr) = cxp(iujvr) should 
be expanded. 

Using power counting and the Z a-expansion, it is easy to 
show that the matrix element in Eq. d55i is of order Za £ 2 , 



which leads to a relative factor of 



a (Z ay 



Z a 



a (Z a) 
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-15 Z 6 



(56) 



for the ratio of the vacuum-induced to the second-order in- 
duced transition probability. The ratio r is very small for 
Z = 1, and it is still only of order 10 -8 for a medium-Z range 
(Z w 20). While it is not negligibly small for heavy ions, the 
transition frequencies increase with increasing Z, and pow- 
erful coherent lasers in this range of frequencies are still un- 
available. 

In order to observe the phenomenon of influence of second- 
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order harmonic generation due to laser-dressed vacuum polar- 
ization on the transition probability, it could be useful explore 
whether transitions exist where the ratio r is more favourable 
with regard to a possible experimental investigation. In this 
case, a very careful estimate of the background would be re- 
quired. For instance, the linewidth of even narrow lasers can 
induce very small, but important transition probabilities far off 
resonance, and these would also have to be taken into account 
in such an analysis. 

D. Application: photon fusion 

Due to the vacuum polarization in the Coulomb field, non- 
linear QED processes such as Delbriick scattering (coherent 
photon scattering in the atomic electric field) and photon split- 
ting (transition of initial photon into two final photons) appear. 
These processes are investigated both theoretically and exper- 
imentally, see Refs. ||T2l [yl Hj, LLJ] • Another closely related 



process is the photon fusion in an atomic field. In this pro- 
cess, two initial photons with the energy u>i and u> 2 transfer 
to the one final photon with the energy lo^ = ui% + u 2 . This 
process has not been observed experimentally because of its 
small probability. 

The existence of the radiation phenomenon has been sug- 
gested by the asymptotic structure of the time-dependent com- 
ponents of the induced scalar potential in Eq. J26l > and the in- 
duced vector potential in Eq. ( I49bl i. Photon fusion takes place 
when two laser photons with the energy to become a single 
photon of frequency 2u>, and this process is otherwise known 
as second-order harmonic generation from the vacuum in the 
presence of a Coulomb field. 

The probability of this process can be calculated as follows. 
First, it is necessary to calculate the electromagnetic field at 
distances r 1/ui, characteristic of radiation. A straight- 
forward calculation leads to the following asymptotics for the 
potential $(r, t) and the vector potential A(r, t): 



$(r,i) 



a(Ze)u 2 
' 720?r m 2 
a(Ze)u) 2 



ml + ~ nl) - 11(£ 2 - £ 2 2 )(n 2 + n\) + 12 i ^n x n 2 



3(£i + e 2 ) - n%) + (ni - n 2 ) (1 - (v ■ ft))] 



exp[2icij(r — t)] 
r(l — (v ■ nj) 



C.C. 



+ 6«66 [{v- n) x («i - rl 2 )} [V ■ (e*i x e 2 )] 



C.C. . 



(57) 



7207rm 2 

- Hfe 2 - £2) + n\) + (fit + n 2 )(l - {v-n))] 

exp[2iw(r — t)} 
?*(1 — (p ■ n)) 

where by C.C. we denote the complex conjugate. Here, we omitted terms proportional to exp[±2i(o; t — ycr)\. These terms lead 
to a energy flux along the propagation direction of the laser field and ensure the fulfillment of the energy conservation condition. 
From these formulas we obtain the corresponding expression for the magnetic field (r> 1 jui), 

i a(Ze)u> 3 
3607T m 2 

-ll(£i-£ 2 )(l-(^))[nx*?] 



H(r,t) 



3(£i + g) [2(fi-u) [v- (nt x n 2 )} + (1 - n))[(rfi - n 2 ) x u]] 



+6i£if 2 [{v- n)(n\ - nj) + (ni - n 2 )(l - {v ■ ft))] [v ■ {e x x e 2 ) 



exp[2iw(r — t)] 



r(l — [y ■ n)) 

The magnetic field is perpendicular to the vector n, as it should be. For the electric field we obtain 

E(r, t) = -V$(r, t) - ^ A(r, t) = [H(r, t) x ft) , 
as should be for the field of radiation. Now we can calculate the differential distribution of the energy flux 



C.C. ; (58) 



(59) 



dS = 



E(r,t) x H(f,t) 



4tt 



r 2 dn 



a(Ze)m 3 

For the total energy flux we have 

2 (a(Ze)uj 



fe + e 2 2 ) 2 + f(c 2 -e 2 2 ) 2 



(n\+^)- l ^t-e 2 ){n\-nl) 



dn 

47T 



s 



1207rm 2 



fe 2 +e 2 2 ) 2 + f(e 2 -i : 



2\2 
2) 



a(Za) 2 uj 2 ( £ 



5400tt 2 \£ 



56 ( &-zr 2 
ai+e 2 



(60) 



(61) 



In view of the suppression factor (£/£o) 4 , this flux is very small, and the observation of this phenomenon demands very 



12 



intense laser sources. Note that the result J6QI > can be also 
obtained with the help of the Heisenberg-Euler effective La- 
grangian (see, e.g., Ref. [XI])- We have checked by explicit 
calculation that our result J60i agrees with the alternative 
derivation using the Heisenberg-Euler effective Lagrangian. 



IV. CONCLUSION 

We have investigated, in detail, the properties of the (one- 
loop) vacuum-polarization tensor in the presence of two non- 
trivial background fields: (i) the atomic, nuclear Coulomb 
field, and (ii) a possibly intense laser field. The potential in- 
duced by the laser-dressed vacuum polarization in a Coulomb 
field is derived in Sec. [H] The approximations employed in 
the calculation concern the frequency of the driving laser ui, 
which we assume to be small as compared to mc 2 /h, where 
in is the mass of the electron [see Eq. 0]. This condition is 
fulfilled by optical wavelengths and for x-rays. Another as- 
sumption concerns the laser field strength £: our calculation 
is valid provided that £ /£q <C 1 [see Eq. dlOi l. where £g is 
the critical Schwinger electric field, but they remain approx- 
imately applicable to situations where, say, £/£q ~ 1/10. 
Consequently, our calculations remain valid for very intense 
laser fields, but they have to be modified when the amplitude 
of the oscillatory electric field is equal or even exceeds the 
Schwinger critical electric field. 

After a derivation of the induced charge density in Sec. Ill Al 
the asymptotics of the laser-induced potential are analyzed in 
Sec. Ill Bl both in a region close to the nucleus (r <C 1/m) 
and in a region which includes the Bohr radius (1/m <C 
r <C 1/uo). The additional laser-induced potential induces en- 
ergy shifts to atomic levels which are analyzed in Sec. Ill CI 
We find that in addition to a familiar 5-state energy shift 
[see Eq. (I32H . which is present even without a laser field, 
there is an additional term which is relevant for P states and 
states with higher angular momenta, which is generated by 
a quadrupole-type contribution to the laser-induced vacuum- 
polarization potential [see Eq. (I35H. The latter contribution 
generates a "polarized Lamb shift" which is akin to similar 
effects predicted in a noncommutative space-time 

The analysis of the charge density and the potential in- 
duced by the laser-dressed vacuum-polarization is only a part 
of the solution of the problem (see Sec. II I It . The laser field 
breaks the rotational symmetry of the pure Coulomb potential 
and induces, in addition, a nonvanishing current density and 
a vector potential (see Sec. MI A> . which gives rise to addi- 
tional fields, which have both time-independent and oscilla- 
tory components. The time-dependent electric and magnetic 
components oscillate at the second harmonic of the incident 
laser frequency and are analyzed in a region which includes 
the Bohr radius of a hydrogenlike system (see Sec. 1111 BV The 
time-dependent electric and magnetic components can lead to 
electric-dipole transitions via second-order harmonic gener- 
ation from the vacuum (see Sec. IHICY In the far-field, the 
induced electric and magnetic fields correspond to the radia- 
tion of real photons, at the second harmonic of the laser (see 
Sec. MIDI . This process is suppressed by a factor (£/£o) 4 



and therefore requires excessively large laser field strengths 
to be experimentally significant. 

Two points are worth mentioning here: (i) all effects de- 
rived in the current investigation are suppressed at least by 
the square of the parameter (£/£q) 2 . This suppression is 
directly related to the fact that the virtual particles in the 
vacuum-polarization loop are effectively "detached" from the 
atomic physics energy scale on which the dominant laser- 
dressed radiative effects are commonly expected 0-0,0,0]. 
This suppression leads to the conclusion that the effects will 
be difficult to observe, experimentally, even at the current 
high-intensity laser facilities. However, with the advent of 
ever more powerful laser facilities which can generate in- 
tense pulses, the effects might become within the reach of ex- 
perimentalists in the not too distant future. On the basis of 
the Heisenberg-Euler Lagrangian, the suppression could have 
been expected, since the additional terms which supplement 
the Maxwell Lagrangian are also suppressed by this factor. 
The Heisenberg-Euler Lagrangian is intimately related to the 
problem studied here (this has been pointed out in Sec. lIIIDV 

One point in particular seems worth a note: The "polar- 
ized Lamb shift" 11611 which is expected in a noncommuta- 
tive space-time, finds a natural counterpart in terms of a laser- 
induced effect which is always present when high-precision 
spectroscopy experiments are performed [see Eq. ( I35H . This 
observation puts a severe restraint on the observability of non- 
commutative space-time effects in atomic spectroscopy. 

In our calculation, we encounter additional finite renormal- 
izations which have to be carried out although the vacuum- 
polarization tensor, which forms the starting point of our cal- 
culation is in principle already renormalized on mass shell 
[see Eqs. Jl 8i and (1421 1 . The reason is that the additional 
external Coulomb field induces spurious effects, which have 
to be eliminated based on the conservation conditions for the 
electric charge and the nonexistence of a magnetic monopole. 
The necessity of carrying out additional finite subtractions is 
akin to the additional subtractions that have to be carried out in 
the calculation of the Wichmann-Kroll higher-order one-loop 
vacuum-polarization corrections in atoms, where, indeed, ad- 
ditional subtractions have to be carried out in each order of the 
Z a expansion. 

Our investigation has been carried out in the more general 
context of quantum field theory under the influence of external 
("unusual") conditions. We have analyzed the vacuum polar- 
ization tensor in the presence of a laser field and a Coulomb 
field, which mediate quantum electrodynamic effects under 
the influence of dynamically changing external conditions. A 
few curious phenomena are encountered: finite renormaliza- 
tions that apply to the induced charge and current density [see 
Eqs. Jl 8i and (I42H . induced charge distributions with an un- 
usual asymptotic behaviour [see Eqs. (I26> . (|42aJ and ( !49bH . 
and the rotational symmetry breaking of the quantum electro- 
dynamic effects [see Eq. ( I35H . The techniques outlined here 
might be useful in other situations where our understanding 
is more limited. In particular, we have shown that general 
wisdom regarding the asymptotic structures of the quantum 
electrodynamic effects has to be modified quite drastically in 
the presence of further external conditions, like a laser field. 
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